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Abstract. Let T be a random triangle in a disk D of radius R (meaning 
f^ . that vertices are independent and uniform in D). We determine the bivariate 

CN i density for two arbitrary sides o, b of T. In particular, we compute that E(a b) = 

(0.837...)-R^, which implies that Var(perimeter) = (0.649...)i?^. No closed-form 
expression for either coefficient is known. The Catalan numbers also arise here. 



PLh 



Let A, B, C denote three independent uniformly distributed points in the disk 
D = m^v) '■ i"^ + ff < R^}- Let T denote the triangle with sides a^h^c opposite 
the vertices A, B, C . We are interested in the perimeter a + 6 + c of triangle T. The 
univariate density /(x) for side a is [Ill2l[3llll|5l|6l[7] 



^ ■ arccos — —^AR^-x"^, < x < 2R 



-kR^ \2RJ 7ri?4 

and 

E(a) = - — R = (0.9054147873672267990407609...)/?, E(a^) = R^. 
457r 

^ '. Clearly 

P : E(perimeter) = 3 E(a) = —R = (2.7162443621016803971222828.. .)i? 

t^^ ■ IStt 

O 

o 



but to compute Var(perimeter) = E(perimeter^) — E(perimeter)^, we will further need 
to consider cross-correlation p between sides. 
The bivariate density f{x,y) for sides a, b is 



><■ f(rv) = S^^^^y^ iix + y<2R, 

^: ■'^'^' \ 'ip{x,y) iix + y>2Randx<2R 

when < y < X (use symmetry otherwise) where 

2xv ( i 

Vi^^y) = ■^\-\/{2R-x-y){x-y){2R + x-y){x + y) + 

o(T^ ^2 ( x^-2Ry + y^ \ , ^^, ( x^ + 2Ry-y^ 
2(n — y) arccos ; : — + 2n arccos 

xy f ft^ + x'-R^\ ft^ + y^-R^\ , 
t arccos arccos at. 



TT^R^ J \ 2tx V 2t|/ 

R-y ^ / \ y 
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^ix,y) = ^y^ tarccos (^ — j arccos (^ j dt. 

It follows by numerical integration that 

E(a6) = (0.8378520652962219016710654.. .)i?^ 

Ii6nc6 

. ,x E(ab) -E(a)E(b) 
p(a,b) = \ ' ^ i ^L^ = 0.1002980835659001715822627..., 

^Var(a) Var(6) 

E(perimeter2) = 3Eia^) + 6E{ab) = (8.0271123917773314100263929.. .)i?^ 

Var(perimeter) = (0.6491289571281667551974101. ..)i?^. 

Exact evaluation of E(o&) remains an open problem. We review derivation of the 
univariate case in section 1, imitating the analysis in [H |9] very closely. (Parry's 
thesis [8] is concerned with triangles in three-dimensional space; it is surprising that 
our two-dimensional analog has not yet been examined.) The bivariate case is covered 
in section 2. An experimental consequence of our work is the formula 

E{aH')='^R' 

which we prove via a different approach in section 3. Finally, in section 4, the Catalan 
numbers from combinatorics appear rather unexpectedly. 

1. Univariate Case 
We omit geometric details, referring to [H |9] instead. The distance t between point C 
and the origin has density 2t/R'^ ioi < t < R. Let f{x\t) be the conditional density 
for distance x between points C and B, given t. We will compute the sought-after 
density f{x) for side a via 

2t 



f{x) = Jf{x\t)—dt. 



There are two subcases. 
1.1. 0<x<R. 

R-x 



, 2x 2t , [ 2x 



R-x 




which corresponds to formula (1.12) in Parry's thesis [8]. The arccos term arises 
since, if the portion of a circle of radius x, center C contained within D has arclength 
26x, then f{x \ t) = {26x)/{itR^); the Law of Cosines gives 6. 
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1.2. R<x<2R. 

f 2x ( e + x^-R^ \ 2i 

x-R ^ ' 

which corresponds to Parry's (1.18). Straightforward integration provides the desired 
result (vahd in both of the preceding regions). 

2. BivARiATE Case 
We omit geometric details, referring to [S] instead. The distance t between point C 
and the origin has density Itj]^ for < t < R. Let f{x,y\t) be the conditional 
density for distance x between points B and C, and distance y between points A and 
C, given t. We will compute the sought-after density f{x,y) for sides a,b via 

f 2t 

f{x,y) = J f{x,y\t)—dt. 



There are six subcases. 

2.1. y < X and < x < R. 

V2x 2y 2t y 2a; ( f + x^ - R^ \ 2y 2t 

R-x ^ ' 

"" 2x fe + x^-R^\ 2y (e + y^-R^\2t ^ 
arccos — — r arccos ^^dt 



7ri?2 \ 2tx / 7ri?2 \ 2ty / R^ 

R~y ^ ^ \ y / 

which corresponds to formula (4.26) in Parry's thesis [8]. 
2.2. R<x <2RandO<y <2R-x. 

y 2x (l+^-R^\ 2y 2t 

x-R ^ ' 

"" 2x (e + x^-R^\ 2y (e + y^-R^\2t ^ 
arccos arccos — -dt 



7ri?2 \ 2tx y 7ri?2 \ 2ty / R^ 

R-y 

which corresponds to Parry's (4.29). Straightforward integration gives ip{x,y) (valid 
in both of the preceding regions). 
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2.3. R < X <2R and 2R-x <y < x. 

^^^'^^ = J^ ^"^^^°^ [ — 2rx — ) ^"^^^"^ [ — 2ry — ) m"^' 

which corresponds to Parry's (4.32). This is, of course, ip{x,y). 

2.4. X <y and < y < R. 

,, , V2x2|/2t, V2x 2|/ /t2 + y2_^2\ 2t , 

^^^'^) = \ WWr^"^'^ J ^;^"^^^°H — ^y — J^ 

R-y ^ y / 

"" 2x (e + x^-R^\ 2y /t^ + y2-R2\ 2t ^ 
arccos — — r arccos --^^dt 



7ri?2 \ 2tx J 7ri?2 \ 2ty I R^ 

R—x 

which corresponds to Parry's (4.35). This is, of course, ip{y,x). 
2.5. R<y <2R amdO <x <2R-y. 

Y2x 2y /t^+jl-R?\ 2t 

tix.y) = / — r- arccos — 7;dt + 

■'^ '^' J R^nR^ V 2ty / R^ 

y-R \ y / 

"" 2x (e + x^-R^\ 2y (e + y^-R^\ 2t , 
arccos — — r arccos -rrrat 



7ri?2 \ 2tx / 7ri?2 \ 2ty / R^ 

R—x 

which corresponds to Parry's (4.38). This is, of course, ip{y,x). 
2.6. R<y <2R and 2R-y <x <y. 

,, , f 2x (e + x^-R^\2y ft^ + y^-R2\2t^^ 

^^^'^) =J^ ^"^^^°^ [ — 2rx — ) ^"^^^°^ [ — 2ry — J ^^' 

which corresponds to Parry's (4.41). This is, of course, ip{y,x). 

3. Characteristic Function 
We follow an approach found in [ini [H]. Let u,v,w denote the squared distances 
between A, B,C and the origin O. Let 99 denote the angle between vectors OA, OB 
and ^ denote the angle between vectors OA, OC. We have 



a = V + w — 2y/vw cos{iIj — ip) 
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b'^ = u + w — 2^/uw cos{^p) , 

(? = u + V — 2^/v/u cos{ip) 

by the Law of Cosines, where -u, w, w are independent uniform on [0, R^] and y9, ip are 
independent uniform on [0,27r]. The characteristic function for {a'^,b'^,c^) is thus 



9{r,s,t) 



1 1 



2tt 2tt 



ij2 r2 f(2 







dip dtp / du dv dw exp ir (v + w — 2^/vw 008(11) — (p)) + 







is (u + w — 2y/uw cos{ip)) + it (u + V — 2y/uv cos{ip) 
It is well-known that 



E(c^) = i $ 



.2\ _ 



dt 



E{b 



.2 ^2 



C^ = -^ 



1 a^^ 



=s=t=0 



dsdt 



r=s=t=0 



and the former becomes 



E{c' 



1 (9 1 1 



27r R2 ^2 



- — -— — / (i(/9 / dudvexp it (u + v — 2y/uv cos{(p) 
i ut R 271 J J J L \ 



1 a 1 ^'^' 

~ ;^"m / / *^-^P (^^ ("" + '^)) "^0 (2tA/m;j du dv 





t=o 



t=o 



where Jo{0) is the zeroth Bessel function of the first kind; hence 

11 r r a 



dudv 



t=0 



E(c') = -^ J J -^exp{tt{u + v))Jo(2tV^\ 



= -^rj / i{u + v)dudv = R^ 
i R* J J 



which is consistent with before. The latter becomes 

''^ "^^ — dip dip I I I du dv dw ex.-^ is (u + w — 2\/uw cos('?/') j + 



E(6'c2 



1 92 1 1 
^dsdtR^A^'^ 







it f -u + f — 2v/m; cos((/9) j 



000 

s=t=0 



d^ 1 
dsdt^ 



i?2 fl2 i?2 



exp 







(is (m + w)) Jo (2sv^-uu'j exp (it (u + v)) Jq \2t^/u/v\ du dv dw 



s=t=0 
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r2 p2 p2 
HUH „2 

~ "r6 / / / a~a7 ^^P ^^'^ ^'^ + ^)) "^^3 (2sy/uw^ exp (it (m + v)) Jq (2ty/mj 



1 111 13^, 



du dv dw 

s=t=0 



— ^ / / I —{u + v){u + w)du dv dw = —R^ 







as was to be shown. The fact that 13/12 — 1 = 1/12 ^ offers the simplest proof we 
know that arbitrary sides of a random triangle in D must be dependent. 

4. Catalan Numbers 
Let i? = 1 for the remainder of our discussion. From 

P (a^ < x) = P (a < v^) = / ( — arccos ( | J - ^y^4-^2J ^^ 



\ TT \'Z I TC 

Q \ \ / 



we obtain that the density for a^ is 



-^^- arccos -^^— V4 - x — ^, < x < 4. 

TT \ 2 / TT / 2^JX 



On the one hand, the characteristic function for a? is 

1 1 

/ / exp {it {u + v)) Jo \2t-\/uv] du dv 



by the preceding section; on the other hand, it is 

f ( ^\fx ( \/x\ X r- \ 1 , 

/ exp(ztx) arccos v4 — x — -^dx 

/ ^ ^V ^ V 2 ; vr^ y 2v/i 

= - [1 - exp(2it) (Jo(2t) - iJii'it))] 

= \[^-Kt)] 

where Ji{9) = —Jq{6). A direct evaluation of the double integral seems to be difficult. 
Boersma |12j, using work of Zernike & Nijboer p^ fl ^ fTS], gave a rapidly-convergent 
series for the inner integral: 

/ exp {itu) Jo {2ty/mj) du = ^^ 1/2 ^^P ("5") ^{-T'T{'in+l)Jn+i/2 f-j ^2™+! (^t^) 
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but this apparently does not help with the outer integral. 

Let Io{0) be the zeroth modified Bessel function of the first kind and Ii{9) = Iq{9). 
We note that the exponential generating function for the Catalan numbers [T6] : 

exp(2«)(/„(20-A(2«))=|^P)r 

is remarkably similar to the expression for h{t). Replacing t by it, we obtain 

°° 1 /2n\ 

hit) = exp(2it) (Jo(2t) - ^M2t)) = ^ -— — (zt)" 

n=o yn-\- I). \n J 

because JoiiO) = Io(0), Ji(i9) = ili(9). Therefore the Catalan numbers are associated 
with the characteristic function for a^. We wonder if a two-dimensional integer array, 
suitably generalizing the Catalan numbers, can be associated with the characteristic 
function for (a^,6^): 

111 



exp (is (u + w)) Jo (2sy/uw] exp {it {u + v)) Jq (2t^/uv] du dv dw. 







Since the bivariate density f{x, y) for (a, b) is much more complicated than the uni- 
variate density f{x) for a, an answer to our question may be a long time coming. 
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